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Abstract 

A cuprate superconductor model based on the analogy with atomic nuclei was shown by lachello 
to have an sn(3) structure. The mean-field approximation Hamiltonian can be written as a linear 
function of the generators of su(3) algebra. Using algebraic method, we derive the eigenvalues of the 
reduced Hamiltonian beyond the subalgebras n(l)0u(2) and so(3) of sm(3) algebra. In particular, 
by considering the coherence between s- and d-wave pairs as perturbation, the effects of coherent term 
upon the energy spectrum are investigated. 
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I. INTRODUCTION 



As is well known, the model of atomic nuclei consisting of s and d nucleon pairs, was 
introduced long ago and has been ver^ successful in describing properties of nuclei 



Recently, it has been suggested by Muller 



that both s- and rf-wave pairing occurs in cuprate 



superconductors with a mixture that varies as a function of the distance r from the surface. 



which stimulated considera 



quantum systems 



considerable mtei 



e interests in studying the novel phenomena of these macroscopic 



A macroscopic model of cuprate superconductors based 



on the analogy with atomic nuclei has been proposed in Ref. jJl, and by developing the OAI 
method Q], lachello pointed out that the basic ingredients of the cuprate superconductor 
model are s- and d-wave pairs in a 2-dimensional Fermi system with a surface, which possesses 
the properties advocated by Muller ^ on the basis of the analysis of recent experiments. The 
analytic solutions of this model based on u{l)(S)u{2) and so(3) subalgebras structures of the 
su{3) algebra have been discussed, which respectively correspond to the case of ignoring the 
interaction strength V2 between s- and d-wave pairs and choosing V2 as some special value, but 
the general solutions going beyond the actual subalgebra-confined framework have not been 
investigated. What about the eigenenergy of the Hamiltonian including the coherence between 
s- and d-wave pairs and the influence of the coherent term on the energy spectra in this system? 
This is our task in this paper. 

Firstly, we study the cuprate superconductors model where there appears to be a compe- 
tition between s- and d-wave pairing, and write its Hamiltonian as the bilinear function of 
the generators of su{3) algebra. Then we use the mean-field approximate (MFA) method 
so that the reduced Hamiltonian of this model can be presented as the linear function of the 
generators of su{3) algebra. Adopting the algebraic method we get the eigenenergies of 
the reduced Hamiltonian by introducing an operator of similar transformation. We also discuss 
the case of considering the coherence between s- and d-wave pairs as perturbation, and analyze 
the effects of coherent term upon energy spectrum subsequently. 

II. THE MODEL AND THE REDUCED HAMILTONIAN 

By mapping from the fermion Hamiltonian, the boson Hamiltonian of cuprate superconduc- 
tors based on the analogy with atomic nuclei was given by lachello jj]: 

H = Ssis^s) + ed{d']^d^ + dtd^) + uo{s~^s^ss) 



+U2{d'Xd'Xd+d+ + d^dtd^d^) + u^i^d'^^d'^d^dJ) 

+Vo{dXd+ + rf+d_)(s+s) + V2{dXdtss + s+s+rf+rf.), (1) 

where Ss and Sd denote the single-boson energies; Uq,U2,u'2,Vq and V2 are the strengths of the 
boson-boson interaction. The creation and annihilation operators of s and d± are denoted by 
s"*", s and c?^, (i± that satisfy the boson commutation relations: 

[s, s+] = 1, [s, d+j = 0, [s, rfj = 0, [du-, dl] = Suu', u = ±. (2) 

An explicit construction of the algebraic generators in terms of creation and annihilation 
operators are: 





= d_dj^. 


, I- 


= dX_d-, 




_ = dX_s, 


U_ 


= s+c?+, 




. = s^d_ 




= dts, 


h 


= ^K- 


- nd_ 


J, 


h 


K 


+ nd_ 


. - 2ns) 



(3) 



where 



Ug = s~^s, Ud = Ud^ + nd_ = + d^d^. (4) 

Based on the relations one can verify that the set {I±, h, I8,U±,V±} satisfies the commu- 
tation relations of Lie algebra su{3): 

[J3, /±] = ±/± [/+, /-] = 2/3 [Is, L] =0 (a = ±, 3) 
[h, U±] = T^Ui [h, U±\ = ±U± [U+, = -(/3 - f/g) 
[h, Vi] = t'2V± [h, Vi] = tV± [V+, y_] = -(/3 + f/g) 
[/±, U±\ = [U±, V±\ = ±1^ [J±, V±] = tU^ 

[I±,U^] = [Ii,V^] = [U^,V^]=0. (5) 
Hence the Hamlitonian can be written in terms of the su{3) generators, i.e., 

H = {2u2 - 1x2)13 + (lio + y + ^ - ^o)^ + [y (-2^0 + 2^2 + ^2 - ''^o) - Ss + ed + uo - U2\Is 

N 1 
+ [y (wo + U2 + 2vo + 2u2) + -{ss + 2ed - uq - 2u2)]N + V2{U+V^ + V+U^). (6) 

Here = + is the total particle number, and [A^, F] = (T = I±,V±,Uj-, I3, Is). 
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The Hamiltonian (jHI) can be reduced to the mean field form by using the standard procedure 
for any bihnear operators such as AB, which can be approximated to the hnear form 

ABc^A<B> + <A>B-<A><B>, (7) 

based on the assumption {A— < A >){B— < i? >) ~ 0, if < A > and < B > are suitable 
expectation values at the ground state. 
Setting the expectation values to be: 

< U+ >= Ai, < >= A*, < V+ >= A2, < V- >= a;, 

< h >= A3, < h >= As, <N>=N, (8) 

and applying the process (jZj) to Hamiltonian (jH)), which leads to the mean- field approximation 
(MFA) of the Hamiltonian ©: 

Hmf = Ho + E', (9) 
Ho = cU+ + dU. + eV+ + fV. + gh + Ws, (10) 

where the coefficients read 

c = f2A2, d = V2A2, e = V2Al, f = V2A1, 
g = 2C/0A3, h = hoAs + hN + h2, 

E' = -t;2(AiA; + AIA2) - goAl - /iqA^ + ciN^ + {e' + hAs + C2)N - hAgN, 

u' 

go = 2u2-U2, /lo = 2wo + M2 + y - 2t;o, 
1 

hi = -(-2mo + 2^2 + U2 - vo), h2 = -Ss + + Mq - "2, 

1 11 

ci = -{uo + U2 + 2vo + 2u2), C2 = --{U0 + 2U2), e' = -{es + 2ed). (11) 

III. THE ENERGY SPECTRUM 

When the Hamiltonian of a model is expressed as a linear function of a Lie algebraic gener- 
ators, there has the algebraic method to diagonalize and obtain the eigenvalues and eigenstates 



ory ^^H]j20|], if Hq is a linear function of the generators of a compact semisimple Lie group. 



of the system. For instance, the algorithm of constructing coherent states |12| and the method 
of algebraic dynamics QlQ] l^E^I were presented. Following the standard Lie algebraic the- 
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it can be transformed into a linear combination of the Cartan operators of the corresponding 
Lie algebra by 

Hi = WHoW-. (12) 

Here W = Dili exp{xiAi) is an element of the group and W~ denotes the inverse of W, in 
which Ai {i = 1, N) is a basis set in Cartan standard form of the semisimple Lie algebra, 
and be set to zero if the corresponding Ai is a Cartan operator. 

In order to give the complete eigenvalues, according to references Q]^^, we choose the 
similar transition operator as 

W = exp{x2iI-)exp{x23U+)exp{x3iV+)exp{x32U^)exp{xisV^)exp{xi2l+), (13) 

and let the coefficients of the non-Cartan operator vanish while substituting Eqs. (jl(jp (|13 p into 
the right-hand side of Eq. lfT^ . then we get a set of algebraic equations of Xij after lengthy 
computation: 

-gxi2 + {d- exi2)xi3 = 

/ + cxi2 -{h+ y)xi3 - exls = (14) 
d - exi2 + {h- + exis - ca;32)x32 = 0, 
and Hamiltonian after the transformation of W becomes diagonal: 

3 3 

Hi = WHqW- = ig + exi3 + cxs2)h + ih + -exi^ - -cx32)h- (15) 

One can see that although operator W is not unitary, the similar transformation (jl2p guar- 
antees that the eigenvalues of Hq equal those of Hi. This is acceptable for we are only concerned 
with the eigenvalues. Supposing the common eigenstates of the Cartan generators and Jg of 
Lie algebra sm(3) are the Fock states | 'nd^,nd_,ns >, i.e., for the commutative set {/a, Jg, A^} 
there exist: 

/3|n.„n._,n. >= l(n,_ - J | >, 

Is I n,„n._,n. >= l(n., + n._ - 2n.) | r..,,n._,n. >, 

N \nd+, nd_ , Us >= {nd+ + nd_ + Us) \ nd+ ,nd_,ns> . (16) 

From Eqs. ()15j)(fTT)|) it follows the eigenvalue of the MFA Hamiltonian Q: 

1 13 3 

Ernf = -{g + exi3 + cx32){nd_ " Ud^) + -{h + -6X13 - -cx32){nd+ + nd_ - 2ns) + E'. (17) 



In the following, we will analyze the solutions of Eq. (jl7j) in detail by comparing with Ref. 



a. 



A. A special case: ^2 = 



Considering the case of V2 = 0, we can easily find that the Hamiltonian Q only consists 
of the su{3) Cartan operators {I3, /g} and the total particle number A^. It is obvious that the 
Eq. ()14|) only has the zero solution for Xij and the eigenvalue reads for V2 = 0: 

E^'^\ns, Ud, I) = EsUs + Band + UQUs^ns - 1) + ^M2^d("'d - 2) + ^u'^nl 

+ {2u2- u'^)!"^ + VortdUs, (18) 

which is the same as "phase (I)" corresponding to m(1)0m(2) algebraic structure in Ref-P]. 
When the system only has s-wave, the eigenvalue is E^^'^'>{ns) = SsUs + UQUging — 1), while for 
p-wave only, the eigenvalue becomes E^^^\nd, I) = SdUd + \u2nd{nd — 2) + ^^^^ + (2^2 — '^'2)1'^ ■ 
(Note that here / = —/a = |(^d+ — which is half of the / defined by lachello \^\-) It 

can be naturally understood that Lie algebra su(3) could reduce to its subalgebra when some 
parameters are specified. 



B. The case of ^2 / 

To discuss the contribution of coherent term in Hamiltonian i.e., the nonlinear term in 
(jni), we will investigate the eigenvalues of the MFA Hamiltonian by considering the coherent 
strength V2 between s- and d-wave pairs as a perturbation. For simplicity, we choose g = ho = 0, 
which means some interaction strengths in Eqs.()lH) satisfy the relations: 

U2 = 2U2, Vq = Uo + U2, (19) 

then the coefficients (lllll in Hamiltonian S reduce to 



0=^2^2, d = V2A2, e = V2Al, f = V2Ai, 

g = 0, h = hiN + h2, 

E' = -V2{AiA; + AIA2) + ciN^ + {e' + /iiAg + C2)N - h^A^N, 

hi = U2- Uo, h2 = -Ss + Ed + Uo- W2, 

1 11 

ci = -{U0 + 2U2), C2 = --(mo + 2M2), e' = -{es + 2ed), (20) 



and the solutions of Eqs. lfT^ read: 



/ 1 

1. Xis = 0, Xi2 = , X32 = — (/i - A), 

c 2c 



/ 1 

2. Xi3 = 0, Xi2 = , X32 = — (/i + A), 

c 2c 
1 d 

3. Xi3 = - — (/l + A), Xi2 = -, X32 = 0, 

2e e 
1 

4. xi3 = - — (/i - A), X12 = -, a;32 = 0, 

2e e 

1 1 

5. a;i3 = — (/i + A), X12 = -, X32 = — (^-A), 

2e e 2e 

6. xi3 = -^(/i - A), .T12 = -, X32 = + A), (21) 

2e e 2e 

where 

A = ^/i2 + 4t;2(AiAl + A2A^). (22) 

Substituting Eq.((2H) into Eq. lfTTjl . we obtain the eigenvalues of reduced Hamiltonian © as: 

/i — 3A /i /i + 3A , 

El = -( — - — )nd+ + - ( — - — )ns + E , 

bib 

^ ,/i + 3A, /i — 3A, ^, 

^2 = - — - — )nd+ + -rid_ - ( — - — )n, + E , 
b 6 b 

h h + 3X h — 3X , 

E3 = -^rid^^ - ( — - — )nd_ - ( — - — )ns + E , 

^ h — 3A, ,/i + 3A, ^, 

E4 = -Ud^ - ( — - — )nd_ - ( — - — )n, + E , 
6 b b 

^ ,h — 3X, ,h — 3X, h ^, 

E5 = - — 7 — - ( — 7 — J'^rf. + -n, + £; , 

D O 3 

„ ,/i + 3A, ,/i + 3A, h . , . 

E^ = -{——)nd^-{——)nd_+-ns + E'. 23 
3 

If coherent coefficient V2 is much smaller than others in the Hamiltonian 0, we regard 
the nonlinear term including V2 as a perturbation one. For f 2 <^ /?. by expanding A to the 
second-order, we get 

A ^ ±h[l + ||(AiAt + A2A;)]. (24) 

Substituting into (PHj) and supposing n^^ = n^. in the following, we obtain 

1. El = -{^)nd^ + - (^)n, + (xi3 = 0, X12 = -f X32 = ^(/i - A)) 
(a) For /i > 0, 



Eia = E^^^ + avl{nd+ - n,) - fiv2, 

= E^'^+avlC^-ns)-f5v2; (25) 
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(b) when h < 




ris) - (3v2 



(26) 



where 




EsUs + SdUd + UQUsius - 1) + U2nd{nd - 1) + {uq + U2)nsnd, (27) 



a = 



and 



(29) 



According to Egs-fpU jlfl^ . /i > 0, or a > 0, corresponds to the case of U2 > Uq 
and Ed > Sg, that is, the energy of s-wave is lower than that of d-wave, and the 
self-interaction strength of d-wave is larger than that of s-wave. On the contrary, 
/i < 0, or a < 0, is ^2 < "^^0 and Sd < Eg, the energy of d-wave is lower than that of 
s-wave, and the self-interaction strength of s-wave is larger than that of d-wave. (3 
denotes the part of the contribution of coherence between s-wave and (i-wave and is 
positive. 

From Eqs.()25 p ()26 |) . one can draw the conclusion that when t'2 = 0, the eigenenergy 
Eia reduces to E^^^ , the eigenenergy of s-wave and (i-wave mixed without the coher- 
ent interaction. Note that when nd+ = nd_ and Eq. (fTTHl holds, the eigenenergy 
(fTHj) reduces to E^^'^ Comparing with Ref.jj|, Ei^ gives another eigenvalue of 

the system, and when V2 = Eu reduces to E^^^\ which is another eigenenergy of 
s-wave and rf-wave mixed without the coherent interaction. Based on Eqs. (j?rj) ()29 p . 
we find that E^^'' > E^^^^ when a > and 2ns < Ud oi a < Q and 2ns > ^d, 
E^^^ < E^^^^ holds for a > and 2ns > Ud or a < and 2ns < nd, i.e., there exists 
the case that E^^^ or E^^^^ is the ground state for different cases. The last two terms 
in Eq. ()25|) and Eq. ()26p are the contributions of the coherent term indicated by f 2 in 
Hamiltonian Q corresponding to eigenenergies Eia and En, respectively, the first 
term including V2 causes split between nd and energy levels, and the last term 
leads to the energy have a total shift, whose direction is determined by the sign of 
V2- Split between nd and and the total shift make the energy spectrum different 
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from the case of V2 = 0, especially the change of the energy gap between s-wave and 
d-wave. 

For h > and ^2 > 0, the energy eigenvalue £"10 has the smaller value when the 
particle number Ug is more than the particle number na, and the coherent strength 

between s-wave and (i-wave is large as possible. For the case of h < and V2 > 0, 
the larger Ug takes and the larger coherence between s-wave and d-wave has, the 
more steady the system is. 

2. E2 = -(^)nrf^ + - (^)n, + E' {x^ = 0, x,2 = -f, 0:32 = ^(/i + A)) 
(a) For h>0, 

h 

E2a = ^"^"^^^ + g (-2nd_^ + nd_ + Ug) - av^na^ + av^Us - /3f2, 

= E^''^-{^ + avl){^-ng)-Pv2. (30) 



(b) When h < 0, 



E2b = E'^^^ + avl{nd+ - Ug) - (5v2, 



= E^'^ + avl{'^-ng)-(5v2. (31) 

Similar to the analysis of case 1, wc can find that for h > and ^2 > 0, the system 
becomes more steady when both na and the coherent strength are larger. And when 
h < 0, the more na and the coherent strength takes, the lower the energy eigenvalue 
has for V2 > 0. 

3. Es - Ina^ - {^)nd_ - (^)n, + E' {x,-, ^-±{h + X), x,2 = f , X32 = 0) 
(a) For h>0, 

Esa = E^^^^ + ^{nd^ - 2nd_ + rig) - avlnd_ + avlug - /?V2, 

= E^''^-{^ + avl){-'^-ng)-Pv2. (32) 



(b) When /i< 0, 



E^b = E^^^ + avl{nd_ - Ug) - pV2, 



Ei^) + avlC-^ - rig) - Pv2. (33) 
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For h > and V2 > 0, the system becomes more steady when both and the 
coherent strength are larger. And when h < 0, the more na and the coherent 
strength takes, the lower the energy eigenvalue has for V2 > 0. 

(a) For h>0, 

= E^'^ + avli"^ - ris) - (3v2. (34) 



(b) When h < 0, 

E4b = E^^^'' + ^{nd+ - '2nd_ + Us) - avlnd_ + avlug - Pv2, 

= E^''^-{^ + avl){-'^-n,)-Pv2. (35) 

When h > 0, the more and the coherent strength takes, the lower the energy 
eigenvalue has for V2 > 0. For h < and V2 > 0, the system becomes more steady 
when both Ug and the coherent strength are larger. 

E5 = -{^)nd+-{^)nd_ + ^n, + E' (^13 = i(/i + A), X12 - f , X32 = ^(/i-A) 

(a) For h>0, 

= + avlud - pV2. (36) 

(b) When h < 0, 

Esb = E^^^^ + ^(-2n<i_^ - 2nd_ + Us) - av^Ud^ + nd_) - Pv2, 
2h h 

= - (_ + avDud + -Us - Pv2. (37) 

Both for /i > and /i < 0, the eigenenergies give new values, and reduce to E^^^^ 
when V2 — 0. When h > and 1^2 > 0, the less Ud has and the larger the coherent 
strength takes, the more steady the system is. For h < and V2 > 0, the more Ug 
has and the larger the coherent strength takes, the more steady the system is. 

Ee^-{!^)nd^-{^)nd_+lns + E' (x^s ^ ^^{h - X) , X12 = f , ^32 = ^(/i + A)) 
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(a) For h>0, 

= -E^^^^ + ^{-2nd+ - 2nd_ + Ug) - avl{nd+ + n^,) - pV2, 

= E^''^ - (y + + ^ris - Pv2. (38) 

(b) When h < 0, 

Eeb = E^^^^ + + avi{nd+ + nd_) - pV2, 

= E^^^^ + avlud - Pv2. (39) 

Both ior h > and h < 0, the eigenenergies give new values, and reduce to E^^^\ 
Either when h > and > or for the case of /; < and > 0, the more has 
and the larger the coherent strength takes, the more steady the system is. 



IV. CONCLUSION 



We have considered a cuprate superconductor model with s- and d-wscve pairs which exhibits 
a dynamical su{3) symmetry. In order to investigate the solutions of this quantum system, we 
have used the mean-field approximation and got the reduced Hamiltonian in a form with linear 
function of the su{3) algebra generators. By introducing the similar transformation operator W 
and eliminating the non-Cartan generators after the similar transformation through parameter 
Xij, we have obtained the eigenenergies of the reduced Hamiltonian beyond the subalgebra- 
confined framework. It should be noted that the order of the operators in W can be chosen 
arbitrarily, but the coefficients Xij are strongly dependent on the order, that is, the operator 
W in this algebraic method is not exclusive. Although any specified order has a solution, a 
properly chosen order can simplify the procedure to get the Xij. The advantage of this algebraic 
method is that the transformation operator W ensures the eigenvalue of Hamiltonian Hi as 
the same as that of Hamiltonian H. In particular, we have analyzed the effects of the coherent 
term, when we consider it as a perturbative one. We get another new eigenvalue E^^^^ for 
the case of s-wave and d-wave mixed without the coherent interaction. Furthermore, the split 
and the whole shift of the energy spectra induced by the coherent term could be understood 
clearly. Much work (such as the numerical calculations) remains to be done in the field of 
searching and seeking more general solutions going beyond the present works, which can help 
us to understand more and deeper about the model of cuprate superconductors based on the 
analogy with atomic nuclei. 
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